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Abstract 

We treat the inverse problem of determining material losses, such as 
cavities, in a conducting body, by performing electrostatic measurements 
at the boundary. We develop a numerical approach, based on variational 
methods, to reconstruct the unknown material loss by a single boundary 
measurement of current and voltage type. 

The method is based on the use of phase-field functions to model the 
material losses and on a perimeter-like penalization to regularize the other- 
wise ill-posed problem. We justify the proposed approach by a convergence 
result, as the error on the measurement goes to zero. 
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1 Introduction 

In many inverse or optimal shape problems arising in the applications, the aim 
is to reconstruct the shape of an object, usually represented by an unknown 
open set, satisfying certain requirements. If we restrict ourselves to a variational 
formulation, for the sake of simplicity, we look for the shape minimizing a given 
functional F among all the admissible shapes. The shape is often modeled as 
a binary function, that is the open set is described through its characteristic 
function. 

Two of the main issues for a satisfactory numerical resolution of this kind 
of problems are the following. First of all, and especially for inverse problems, 
the problem may be ill-posed, that is stability is missing or, in other words, 
F is not continuous. Second, numerically handling shapes or sets is not an 
easy task from the implementation point of view. The first issue is usually 
tackled by a regularization method, namely by adding to the functional a term 
penalizing the binary function with respect to some -BT^-related norm. For most 
applications, this should be enough for ensuring a regularization without being 



restricting on the class of admissible unknowns. Often the -BF-related norm is 
simply a perimeter-like penalization. About the handling of shapes or sets in 
computations, in many cases this is performed by associating to the open set 
a smooth function describing it. For example, one way of doing it is to replace 
the characteristic function of an open set D with a smooth function, referred to 
as a phase-field function, which is close to outside D, close to 1 inside D, and 
has a quick transition from to 1 across the boundary of D. Another way is the 
so-called level-set method, where D is identified with the sublevel set {tp < 0} 
of a smooth function ip. 

We are interested in using perimeter-like regularizations and phase-field func- 
tions for solving inverse or optimal shape problems, in particular those that are 
not well-posed. We aim to prove in a rigorous way that this kind of approach 
provides a good approximation of the original problem, allowing us at the same 
time to tame the ill-posedness and to have a formulation amenable to be easily 
implemented. A cornerstone of this method is the approximation, in the sense of 
T-convergence, of the perimeter functional by functionals defined on phase-field 
functions, due to Modica and Mortola, [TB]. Since [TS] such a result has found 
innumerable applications. In fact, whenever the functional F is continuous in a 
suitable way, the invariance of T-convergence by continuous perturbations per- 
mits to obtain an analogous T-convergence result if we add to F the perimeter 
penalization. Whenever the problem is ill-posed, that is F is not continuous, 
a corresponding convergence result is not straightforward any more. Since we 
believe that the method is valuable also in the ill-posed case, it would be im- 
portant to justify it in a rigorous way, in general through a convergence result 
inspired by T-convergence techniques, for various interesting applications. 

In this paper we perform such an analysis for the following inverse problem, 
arising from non-destructive evaluation. We aim to determine perfectly insu- 
lating defects in a homogeneous and isotropic conducting body by performing 
electrostatic measurements of voltage and current type at the boundary. The 
conducting body is contained in f2, a bounded domain of M. N , N > 2. The de- 
fects may have different geometrical properties, for instance we may have at the 
same time cracks (either interior or surface breaking), or material losses (either 
interior, that is cavities, or at the boundary). We denote with K the union of 
the boundaries of these defects, whereas 7 is a part of the boundary of f2 which 
is accessible, known and disjoint from K. If a current density / € L 2 (7), with 
zero mean, is applied on 7, then the electrostatic potential u — u(f,K) is the 
solution to the following Neumann boundary value problem 



We call Gk the connected component of £l\K which is reachable from 7 and we 
say that a defect is a material loss if Gk is equal to the interior of its closure, 
that is if no crack-type defect is present. 

The value of u, that is the voltage, may be measured on another part of 
the boundary of ft, say 7, which we assume to be accessible, known, disjoint 
from K and belonging to BGk- We call g such a measurement, that is g = u| 7 . 
For simplicity, we may also assume that 7 coincides with 7. If the defect K is 
unknown, we aim to recover its shape and location, that is Gk, by prescribing 



(1.1) 




Au = in n\K, 
Vw ■ v = f on 7, 
Vw • v = on d(n\K)\j. 
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one or more current densities / and measuring the corresponding voltage on 7, 
g = u\-y, where u solves (jl.ip . In mathematical words, we are given one or more 
pairs of Cauchy data (g, /) on a known part of the boundary and we aim to 
reconstruct the domain of validity of the elliptic equation. 

Here we are interested in the reconstruction only of material losses, that is 
cavities or material losses at the boundary, and for simplicity we refer to it as 
the inverse cavity problem. It is well-known that, in every dimension, a single 
boundary measurement is enough to reconstruct a material loss, thus providing 
uniqueness for the inverse problem, see for instance [211 for a proof with minimal 
regularity assumptions on the unknown material loss. Stability results have been 
proved in [1] for the three dimensional case and in [2] for the planar case, where 
also the instability character of the problem has been explicitly shown. 

We notice that u, the electrostatic potential solution to (11.11) . is constant 
on any connected component of Q\K different from Gk ■ The key observation 
is that its jump set in f2 is essentially contained in K . The uniqueness result 
recalled before actually allows us to say more, in fact the jump set of u uniquely 
determines Gk, that is the unknown material loss. Therefore we are interested in 
the reconstruction of the electrostatic potential u and especially of its (unknown) 
jump set. This suggests the possibility to set up a reconstruction procedure by 
solving a free-discontinuity problem related to the function u. 

The main difficulty for the reconstruction is due to the ill-posedness of the 
problem. In fact, since they are measured, the Cauchy data that are available 
are not exact. Since the problem is severely ill-posed, such an error on the 
measurements may lead to a much greater error on the reconstructed defect. 
Furthermore, the inverse problem is nonlinear. In fact, even if the direct prob- 
lem (jl.ip is linear, the dependence of the electrostatic potential u, and of its 
values on 7, from the defect K is nonlinear. Finally, from a numerical point of 
view, the fact that the unknown is a set, namely Gk, introduces an additional 
complication for the implementation. 

We propose a variational method to tackle at the same time these difficulties. 
The idea is to use a perimeter-like penalization to regularize the problem and 
to replace the unknown set Gk with its characteristic function and, in turn, 
with a phase-field function, to obtain a formulation that may be implemented 
numerically. Namely, the regularization we propose is related to the so-called 
Modica-Mortola functional, an approximation of the perimeter when phase-field 
functions are used. We might construct a family of functionals, depending on 
the noise level on the measurements e, to be minimized with respect to the 
variable u (the reconstructed potential) and the phase-field variable v. However, 
to simplify the implementation we would rather have a functional depending on 
the phase-field variable v only. Thus, we take u depending on v, u — u(v), as 
a solution to an almost degenerate elliptic problem whose coefficient is given 
by a slight modification of v, depending on s. In other words, we replace the 
direct cavity problem with an elliptic problem in where the coefficient of the 
equation is close to 1 in Gk, close to outside Gk, with a quick transition across 
the boundary Gk- The method consists then of minimizing the so-obtained 
functionals, depending on e, with respect to the phase-field variable v only. 
We remark that the reconstructed material loss may be simply computed by a 
suitable thresholding of the minimizing phase-field and that an approximation 
of the looked-for electrostatic potential is given by it = u(v) where v is the 
minimizing phase-field. 
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The main result of the paper, Theorem 14.21 is that the corresponding min- 
imizers v £ converge, as e — > + , to the characteristic function of Gk, thus 
identifying the looked- for material loss, and that u £ — u(v £ ) converge to the 
looked-for potential u. Such a convergence result, whose proofs is obtained by 
techniques borrowed by T-convergence, provides a rigorous justification of the 
method. About the material loss to be reconstructed, this is assumed to satisfy 
a Lipschitz type regularity. We finally remark that the method makes use of 
a single measurement and that is enough to reconstruct the whole unknown 
material loss K. 

If we instead allow the unknown defect not to be a material loss, that is it 
may include crack-type defects, for simplicity we refer to this problem as the 
inverse crack problem. About uniqueness, stability and reconstruction results on 
the inverse crack problem, we refer to [7] and the references therein. The main 
difference between the two cases is that for the determination of cracks one mea- 
surement is not enough, however, at least in the planar case, two suitably chosen 
measurements are sufficient. In [22l [23] a corresponding variational approach for 
the inverse crack problem has been developed. Again such an approach makes 
use of a penalization on the (N — l)-dimensional measure of the defects and of 
phase-field functions. Namely, it was constructed a family of functionals, again 
depending on the noise level on the measurements e, to be minimized with re- 
spect to the variable u (the reconstructed potential) and the phase-field variable 
v. Instead of the perimeter functional and the Modica-Mortola functional, the 
Mumford-Shah functional [TB] and its approximation, in the sense of T conver- 
gence, due to Ambrosio and Tortorelli [H[S] were used, respectively. Also in this 
case a convergence result guaranteed a justification of the method. In Section [5] 
we recall the results obtained in [2 3) for the inverse crack problem and we com- 
pare with those obtained here for the inverse cavity problem. The main difficulty 
in the implementation of the method of [53] is that the functional to be mini- 
mized depends on two variables, the variable u, which should approximate the 
electrostatic potential, and the variable v, which is the phase-field variable that 
should approximate the jump set of the potential and hence the defect. It would 
be desirable to formulate the problem depending on one variable only, for in- 
stance only on the phase-field variable. Unfortunately such a formulation, which 
is proved here for the material loss case, may not be feasible. In fact, Section [5] 
is devoted to show that the result in [33J is essentially optimal, through several 
counterexamples. Moreover, more regularity is needed for the unknown defects 
of crack-type, namely a regularity assumptions of C 1 type, instead of Lipschitz, 
have to be imposed. Thus, we show that restricting ourselves to the reconstruc- 
tion of material losses allows us to gain the following advantages. First, we may 
lower the a priori assumptions on the unknown defect to ones which are more 
suited for applications. More importantly, we obtain and justify a formulation 
which looks more natural and quite simpler to be implemented. 

We finally wish to mention that a numerical implementation, based on the 
results of this paper and on those of [53J , may be found in [20J . The corresponding 
numerical experiments show the validity of these methods also from a practical 
point of view. 

The plan of the paper is the following. After a preliminaries section, Sec- 
tion [2] we describe the setting of the direct and inverse problem in Section [3] 
We treat the material loss case in in Section 21 where there is the main result 
of the paper, Theorem 14.21 In Section [5] we recall the results for the inverse 
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crack problem proved in J23 and we compare the crack and material loss cases 
and discuss their differences. In particular we show the optimality of the result 
of [53] . Finally, in Section [5] we deal with the differentiability of the functionals 
involved. Such differentiability is crucial for developing the algorithm used in 

can- 
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2 Preliminaries 

Throughout the paper the integer N > 2 will denote the space dimension. We 
remark that we shall sometimes drop the dependence of any constant upon 
N, the space dimension. For every x £ l w , we shall set x = (x',xn), where 
x' £ K* -1 and xn & R, and, for any r > 0, we shall denote by B r {x) the open 
ball in 1 N centred at x of radius r. Usually we shall write B r instead of B r (0). 
For any subset E C and any r > 0, we denote B r (E) = {J xGE B r (x). 

For any non-negative integer fc we denote by T-L k the fc-dimensional Hausdorff 
measure. For Borel subsets of M. N the iV-dimensional Hausdorff measure coin- 
cides with C N , the iV-dimensional Lebesgue measure. Furthermore, if 7 c l w 
is a smooth manifold of dimension fc, then T-L k restricted to 7 coincides with its 
fc-dimensional surface measure. For any Borel E C R N we let \E\ — C N {E). 

We recall that a bounded open set 51 C is said to have a Lipschitz 
boundary if for every x £ <951 there exist a Lipschitz function tp : R^ -1 — > K 
and a positive constant r such that for any y £ B r (x) we have, up to a rigid 
transformation, 

y G 51 if and only if yjq < ip(y'). 

We observe that in this case the boundary of 51 has finite (N — l)-dimensional 
Hausdorff measure, that is W. N ~ 1 (dn) < +00. 

We say that a function <p : A — > B, A and B being metric spaces, is bi- 
Lipschitz if it is injective and ip and (p^ 1 : ip(A) — > A are both Lipschitz func- 
tions. If both the Lipschitz constants of <p and c/? -1 are bounded by L > 1, then 
we say that p is bi-Lipschitz with constant L. 

We recall some basic notation and properties of functions of bounded varia- 
tion and sets of finite perimeter. For a more comprehensive treatment of these 
subjects see, for instance, [3J [TTJ [T3] . 

Given a bounded open set 51 C M. N , we denote by BV(ft) the Banach space 
of functions of bounded variation. We recall that u £ BV(tt) if and only if 
u £ L 1 (il) and its distributional derivative Du is a bounded vector measure. We 
endow BV(ft) with the standard norm as follows. Given u £ BV(Cl), we denote 
by \Du\ the total variation of its distributional derivative and we set ||u||sv(f2) = 
IMIli(o) + \Du\(Cl). We recall that whenever u £ W 1 ' 1 ^), then u £ BV(Q) and 
\Du\(fl) = j n \Vu\, therefore ||u||sv(n) = IMUi(fi) + I|Vm||li(o) = \\u\\ w i,i(n). 

We say that a sequence of BV(Q) functions {uf l }'^L 1 weakly* converges in 
BV(il) to u £ BV(ft) if and only if Uh converges to u in L 1 (51) and Duh weakly* 
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converges to Du in Q, that is 

(2.1) lim / vdDuh = / vdDu for any v G Co(f2). 

h Jn Jn 

By Proposition 3.13 in [3], we have that if a sequence of BT/(f2) functions 
{ u h}hL\ is bounded in BV{VL) and converges to u in L 1 ^), then u G BV(il) 
and it/j converges to u weakly* in BV(fl). 

Let Obea bounded open set with Lipschitz boundary. A sequence of BV(Q) 
functions {uh}f? =1 such that sup h \\uh\\BV(n) < +°° admits a subsequence con- 
verging weakly* in BV(fl) to a function u G -BV(f2), see for instance Theo- 
rem 3.23 in [3]. As a corollary, we infer that for any C > the set {u G .BV^fi) : 
ll M ll w(n) < C} is a compact subset of L x (r2). 

Let E be a bounded Borel set contained in Br C R N . We shall denote by xe 
its characteristic function. We notice that E is compactly contained in Br + i, 
which we shall denote by E <g Br + \. We say that -E is a set of finite perimeter 
if belongs to BV(Br + i) and we call the number P(E) — \Dxe\{Br+i) its 
perimeter. 

Let us further remark that the intersection of two sets of finite perimeter 
is still a set of finite perimeter. Moreover, whenever E is open and H N ~ 1 (dE) 
is finite, then E is a set of finite perimeter, see for instance [TTJ Section 5.11, 
Theorem 1]. Therefore a bounded open set f2 with Lipschitz boundary is a set 
of finite perimeter and its perimeter P(£l) coincides with H N ^ 1 (dH,). 

For any bounded open set Q, we define the following perimeter functional 
P : L 1 ^) -> [0, +oo] such that 

!c\Du\(Sl) if u G BV{VL) and u G {0, 1} a.e., 
+ oo otherwise, 

where c is a positive constant to be chosen later. We observe that P(u) = cP(E) 
if u = xe and E is a set of finite perimeter compactly contained in Q. 

We denote by SBV(£l) the space of special functions of bounded variation. 
For any u £ SBV(ft), the density of the absolutely continuous part of Du 
with respect to C N will be denoted by Vu, the approximate gradient of u. The 
singular part, with respect to C N , of Du is concentrated on J(u), J(u) being 
the approximate discontinuity set (or jump set) of u in VL. We further say that a 
function u G GS , SV(f2), the space of generalized functions of bounded variation, 
if it G L 1 ^) and for any T > its truncation u T = (-T) V(TAu) G SW(ft). 
Let us recall that the approximate gradient Vw of u € GS-EiV^ri) is defined 
almost everywhere and coincides with Vmt almost everywhere on {u — ut}, 
and that J(u) = Ut>o J{ u t)- 

The special functions of bounded variation have important compactness and 
semicontinuity properties, see for instance [3] Theorem 4.7 and Theorem 4.8]. 

We remark that if u G BV(SY) and u € {0, 1} almost everywhere in ft, then 
u G SW(fJ) and P(ti) = c|L>w|(0) = cH A T ~ l (J (u)) . 

Let us define the so-called Mumford-Shah functional, introduced in |18] in 
the context of image segmentation. Let us fix positive constants b and c. Let 
MS : L 1 (0) -> [0, +oo] be given by 

(2.3) MS(u) =b [ \\7u\ 2 +cH N - 1 (J(u)) 'due GSBV(Q), 
Jn 



6 



whereas A4S(u) = +00 otherwise. 

Let us introduce the following T-convergence results concerning the approxi- 
mation of the perimeter functional and the Mumford-Shah functional by phase- 
field functionals. For the definition and properties of T-convergence we refer to 
[8]. The perimeter approximation is due to Modica and Mortola, [16], whereas 
the Mumford-Shah functional approximation is due to Ambrosio and Tortorelli, 
[H [5] . We shall follow the notation and proofs contained in [5] . 

Throughout the paper, for any p, 1 < p < +00, we shall denote its conjugate 
exponent by p', that is p _1 + (p')^ 1 = 1- Let W : R — > [0, +00) be a continuous 
function such that W(t) = if and only if t G {0, 1}. Let c w = f ^W{t)dt. 
In the definition of the perimeter functional we pick c = 2cw- For instance, we 
may choose W(t) = 9t 2 (t — l) 2 for any t £ R, whence c = 2cw = 1. 

The following approximation result is due to Modica and Mortola, [IB], see 
also [6] Theorem 4.13]. 

Theorem 2.1. Let fi C R be a bounded open set with Lipschitz boundary. 
For any rj > we define the functional P v : L (O) — > [0, +00] as follows 

!- I W(v)+r] ( |Vw| 2 ifve W 1 ' 2 ^, [0, 1]), 
V Jn Jn 
+ 00 otherwise. 

Then we have that, with respect to the metric of L 1 ^), P, q T-converges to 
P as 7} ->• 0+ . 

Here W 1 ' 2 ^, [0, 1]) = {v £ W 12 {Q.) : < v < 1 a.e. in O}. We note that 
the result does not change if in the definition of P v we omit the constraint 

< v < 1 a.e. in fi. 

Also the following result, due to Modica, [15] , will be useful. 

Proposition 2.2. For any C > and any r\ > 0, let us define 

A c = {v £ L x (0) : < v < 1 a.e. and P v (v) < C}. 
Then Ac is precompact in L 1 (f2). 

Remark 2.3. With the same proof, we can show the following. Let us consider 
any family {w^}o<i ) <i )0 such that, for some positive constant C and for any 77, 
< 77 < rjo, we have < v v < 1 almost everywhere and P^(u^) < C. Then 
{ v v}o<n<no i s precompact in 

Let us fix q, 1 < q < +00. Let V : M — > [0, +00) be a continuous function such 
that V(t) = if and only if t = 1 and let c y = >/y(t)dt. Let V : K -> M be a 
continuous non-decreasing function such that ip(0) — 0, ^(1) = 1 and tp(t) > 
if t > 0. For any 77 > 0, let us fix o v = o v (q) > such that lim^^ + °r]M q = 0. 
Finally, we define ipj, = (1 — o v )ip + o n . Provided o n < 1, we have that ip^ is a 
continuous, non-decreasing function such that ip v {0) = o v and ip v 0-) = 1- 

For instance, we may choose V(t) = (t — l) 2 /4 for any t £ R, whence 
4cy = 1. About tp, we may take ip(t) = f , 7 > 0, if t > 0, while ip(i) = if 
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t < Alternatively, we may choose ip(t) = if t < 0, ^(t) = — 2t 3 + 'it 2 for any 
t £ [0,1], and ^(t) = 1 for any £ > 1. We may finally take o ri (q) = r/ q . 

Then, for any n > 0, we define the following functional AT q n : i 1 (r2) x 
L x {Vl) [0,+oo] by 

(2.5) .A7l(u,?;) = & / ^»|Vw|« + - / V(v)+ V f \Wv\ 2 
Jn V Jn Jn 

if u G W 1,q (Q) and v £ W 1<2 {n, [0, 1]), 

whereas AT^(u, v) — +00 otherwise. We shall refer to AT q as the Ambrosio- 
Tortorelli functional. 

Let us define the following variant of the Mumford-Shah functional. The 
main difference is that we allow the exponent q to be different for 2, requiring 
only that 1 < q < +00. For reasons which will appear evident soon, we also 
add a formal variable v and we pick c = 4cy. We define the functional A4S q : 
x L^ft) [0,+co] by 



(2.6) MS q (u,v)=b \Vu\ q +4c v H N ~ 1 (J(u)) 



N-l, 

v u| - T 'i^v' n. 1 

n 

if it G GSW(ft) and u = 1 a.e. in ft, 

whereas A^5 9 (u,w) = +00 otherwise. 

The Ambrosio-Tortorelli functional approximates the Mumford-Shah func- 
tional, in the sense of L-convergence. Such an important approximation result 
is due to Ambrosio and Tortorelli, [HE], see also Ej. 

Theorem 2.4. With respect to the metric o/L 1 (J7) x L 1 (J7), we have that, as 
f] -> 0+, AT 9 , r -converges to MS q . 

Let us review some regularity results which will be needed in the sequel. 
Most of these results are a consequence of a theorem by Meyers, [17], see also 
(12j . and of standard regularity estimates, and we shall omit the proofs. Let 
ft C be a bounded open set with Lipschitz boundary. Let A = A(x), x £ ft, 
be an N x N matrix whose entries are measurable and such that, for some 
< A < 1, we have 

, _s A[x)i ■ £_ > A|£| 2 for any £ £ R. N and for a.e. x G ft, 

(2 - 7) < A- 1 . 

We remark that for any matrix A, by || A\\ we denote the norm of the matrix as 
a linear operator. 

Let / £ I. s (<9ft), with s > 1 if N = 2 or s > 2(N - 1)/N if N > 3, 
be such that J 9Q f = and let F G L^ft,]^), with p > 2. Let us denote 
W*' 2 (ft) = {116 W 1 ' 2 ^) : J n u = 0}. Then, there exists a unique u G Wl' 2 {ti) 
such that 

(2.8) / AVu-Vv= I F-Wv+ I fv for any v £ W 1,2 (ft). 

Jn Jn Jon 

This is the weak formulation of 

div(AVu) = div(F) in ft 
AVu ■ v = / on 9ft. 

The following regularity result holds true. 
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Proposition 2.5. Under the previous assumptions, the following regularity 
properties hold. 

First of all, we have, for a constant Co depending on N, X, p, s and f2 only, 
(2.9) \\u\\ wl ,2 {Q) < Co (\\F\\ LP{ n) + H/lkW • 

If p > N and s > N — 1, there exist a constant C\ > such that 
(2-10) IMU~ (n) < Ci {\\F\\ LP{n) + \\f\\Ls { en)) ■ 

Here C\ depends on N , X, p, s and fi only. 

There exists a constant Q > 2, depending on N , X and Q only (Q — > 2 if 
X — > + ), such that if p satisfies 2 < p < Q and s > (N — l)p/N, then 

(2.H) I|V U |Up ( o) < C 2 {\\F\\ LP{n) + \\f\\ L . {B a)) ■ 

Here Ci depends on N, X, p, s and f2 only. 

We conclude that if s > N — 1, there exists a constant q(X) > 2, depending 
on N, X, s and f2 only, such that for any p, 2 < p < q(X), we have 

(2-12) HVull^n) < C 3 (\\F\\ LPm + \\f\\ L . (Bn) ) , 

in particular, if p = q(X) 

(2.13) \\Vu\\ LgW{n) < d (||F|| ilW(n) + ||/|U. (a n)) • 

Here C3 depends on N , X, p, s and only, whereas C4 depends on N, X, s and 
fl only. 

Remark 2.6. Let us observe that Q and q(X) converges to 2 as A — > + , whereas 
all the constants C0-C4 might tend to +00 as A — > + . Let us also remark that 
the same kinds of estimates hold true if we replace Wl' 2 {VL) with, for instance, 

W 1 /^) = ju e W 1>2 (0) : J m = o| 

where E is a Borel subset of dtt with non-empty interior, clearly with respect 
to the induced topology of dfi. In this case, the constants C0-C4 might depend 
on E as well. 

We conclude this section with the following lemma, in which we state a 
Caccioppoli inequality. 

Lemma 2.7. Let us assume that A = A(x), x G P>2R, is a symmetric N x 
N matrix whose entries belong to L°°(B2r). We also assume that, for some 
constants < A < A, we have 

A|£| 2 < A(x)£ ■ C < A|^| 2 for every £ G R N and for a.e. x e B 2R . 

Let w G L°°{B2r) be a weight satisfying < e < w < 1 almost everywhere in 
B2R- If u solves in a weak sense 

div(wA\7u) =0 in B 2R , 

then 

(2.14) f w\Vu\ 2 < -2 / wu 2 

Jb r H J B2R 

where C depends on X and A only. 
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Proof. In order to prove (|2.14j) it is enough to take a cutoff function x such 
that x S G£ d (B 2 r), < x < 1 on B 2 r, and x = 1 on Br. We may also assume 
that for an absolute constant C we have |V%| < C/R on B 2 r. Then we use 
the test function ux 2 in the equation and with simple computations we obtain 

dm. □ 



3 The direct problem and setting of the inverse 
problem 

Let 0, Sli and Oi be three bounded domains contained in R , N > 2, with 
Lipschitz boundaries such that Oi C Oi C O and the following properties are 

satisfied. First, O\0x is not empty and dist(Oi, 90 1 HO) > 0. Then, there exist 
7 and 7, closed subsets of 90, which are contained in the interior of 90 n 90i 
and whose interiors, with respect to the induced topology of 90, are not empty. 

We assume that O, Oi, Oi, 7 and 7 are fixed throughout the paper. 

Let Kq be an admissible defect, that is Kq is a non-empty compact set 
contained in O such that dist(iQ)i Oi) > 0. We denote with Gk the connected 
component of £1\Kq such that Oi C Gk - We observe that 7 U 7 C 9Gav We 
remark that if Ko C 30 then Ga' = O, that is no defect is present in the 
conductor. 

We say that an admissible defect Kq is a material loss defect, or material loss 
for short, if Gk is equal to the interior of its own closure (that is no crack-type 
defect is allowed). 

Remark 3.1. If the defect Ko to be reconstructed is compactly contained in 

and 90 is connected, then we may take 7 and 7 equal to 90. Furthermore, 
if 90 is regular enough and we a priori know that dist(iio, 90) > 5, for some 
5 > small enough, then we may choose Oi = {x € O : dist(a;, 90) < 6/2} and 

01 = {x G O : dist(x, 90) < 35/ A}. 

Let us fix a number s, s > N — 1, which shall be kept fixed thorughout 
the paper. Let us prescribe /o G L s (dil) such that j da fo — 0, /o ^ and 
supp(/ ) C 7. 

Let the electrostatic potential uq — u(fo,Ko) be the weak solution to the 
following Neumann boundary value problem 

!Au = in fl\K , 

Vtio ■ v = f on 7, 
Vu • v = on 9(O\if )\7, 

with the normalization conditions 

(3.2) [ u = 0, 

and 

(3.3) uq = almost everywhere in O\G^- . 

Let us recall that our measured additional information is go G L 2 {^) where 
go = uq| 7 . By (EH), we have J g = 0. 
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We observe that there exists a unique solution uo to (I3.ip - p.2p ~ p. 31) and 
that it satisfies the following regularity properties, see [22J for further details. 

There exists a constant C\ > 0, depending on s, f2, f2i, Qi, 7 and 7 only, 
such that 

(3-4) ||V«o|U- (n \jfo) < Ci||/olU-(7). 

(3.5) ||«o||i-o(n) < Ci||/o|U-( 7 )- 

The estimate p.5p guarantees that uq belongs to W 1,2 (0\Kq). Furthermore, 
under the additional assumption that r H N ~ 1 (K ) < +00, or equivalently that 
H N ~ 1 (dGK ) < +00, we have that uo belongs to SBV(£l), its approximate 
discontinuity set J(u ) satisfies H N ~ 1 { J(uo)\dGK ) = anc ~, finally, Vwo, the 
weak derivative of uq in Q\Kq, coincides almost everywhere in fi with the ap- 
proximate gradient of uq, see for instance [5J Proposition 4.4]. 

For any r, 1 < r < +00, and any Borel set E C dil whose interior, in the 
induced topology, is not empty, we define 

w 1 /^) = L w^ r (n) : / u = 



We observe that, by a generalized Poincare inequality, on W E ' r (Ct) the usual 
W 1:r (^l) norm and the norm ||«||vir L ' r (n) = ||Vu||i,r(n) are equivalent. Therefore, 

we shall set this second one as the natural norm of W'/^i). 

Let us consider a weight w in Q satisfying the following properties. We 
assume that w G L°°(SY) and that w > e almost everywhere in 17, for some 
e > 0. 

For any such weight w, and any u%, U2 G W 1,2 (f2), we define the bilinear 
form 

("Ul,M 2 )to= / wVlii ■ V« 2 



and we denote the seminorm 

1/2 



\ux\ w = {ui,ui)l/ 2 = yj w\Vui\ 2 
We denote, for any u\ G IF 1,2 (17), 

HwilU = (||wi||i2(fi) + 



1/2 



We have that || • is an equivalent norm for W 1,2 (tt), and (•,•}«, is a scalar 
product on W E ' (17) whose corresponding norm, | • \ w , is an equivalent norm for 
Wg' (O), for any set E as before. 

For any such weight w and any / G L s (dfl) such that J gn f — and 
supp(/) C 7, let u = u(w) be the solution to the following Neumann type 
boundary value problem 

div(wVu) =0 in 17 
(3.6) { wVu ■ v = f on dft 
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The weak formulation of (|3.6p is the following. We look for a function u <E 
W*' 2 (fi) such that 

/ aiVu ■ Vui = / fui for any u\ € lF 1,2 (r2). 
Jn Jj 

Obviously we have existence and uniqueness of such a solution. Furthermore, 
the following regularity result holds for u. 

Proposition 3.2. Under the previous notation and assumptions, let u solve 
for some weight w. We assume that 

IIHU-(fi) <a 



w{x) — 1 for a.e. x € fij.. 

Then there exists a constant C2, depending on s, 57, fii, fii, 7, 7 and A 
only, such that 

(3-7) \u\ w < C 2 ||/|U^) 

(3.8) HU~(n) < C 2 \\fh»(y). 

We notice that the constant Ci does not depend on w or on e. 



Proof. We sketch the proof of this proposition. Inequality (|3.7p follows from an 
application of Poincare inequality in Ox- The L°° bound (|3.8p is a consequence 
of the maximum principle and may be proved following the same arguments 
used to prove (|3.5I) . see [35] for details. □ 

Let us fix the notation for our inverse problem. Let Kq be the unknown 
defect, which for the time being we assume to be just an admissible defect. 

We assume that /o belongs to L s {d£l) and satisfies supp(/o) C 7 and 
fdn f° = 0- We recall that s is a fixed constant such that s > N — 1. 

The unknown electrostatic potential is u — u(K ,f ), solution to (I3.1[) - 
and the additional measured data is go = wq| 7 . We observe that 



.go e £ 2 (7) and J y g = 0. 

Let us fix a noise level e, < e < 1, then the noisy Cauchy data are given by 
f e and g £ . Here f e belongs to L s (dil) and satisfies supp(/ e ) C 7 and J dn f e = 0, 
whereas g e belongs to L 2 (j) and satisfies J^g £ = 0. We assume that 

(3-9) ||/o - f e \\L*(y) < £ and \\g a - g e \\ L 2 (l) < e. 

For any < e < 1, let 77 = i](e) > and a E > be such that lim e _ s . + 77(e) = 
and lim e ^ + a e = 0. Further assumptions on 77(e) and a e will be imposed later. 

Let us fix a constant ci, < c\ < 1. We recall that tp : M — > K is a continuous, 
non-decreasing function such that ^(0) = 0, t/>(1) = 1, and ^p{t) > if t > 0. 
In particular ip( c i) > 0- Provided < 1/2, we have that ip^ is a continuous, 
non-decreasing function such that 7/^(0) — o n and tpnO-) = 1- Furthermore, 
^(ci) > ^(ci)/2 > 0. 

In the sequel we shall always assume that 

< o n < 1/2 for any 77 > 0. 
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Without loss of generality, we also assume that ip, W, and V are bounded all 
over R, for instance by a constant A. For any r\ > 0, again without loss of 
generality, we assume that ip is such that ip v > o v /2 all over R. 

To any function v £ we associate the function v = 1 — v. We observe 

that, provided < v < 1 almost everywhere in f2, we also have < v < 1 almost 
everywhere in f2. 

For any r\ > and for any v £ L 1 (f2), let = = ip v (v), where 

v = 1 — v. We observe that is such that Hu^l^Dcm) < A + 1/2 and > o v /2 
almost everywhere in Q. Therefore we define, for any v £ L 1 ^) and any e, 
< e < 1, the function u £ £ Wy 2 (f2) where u e = u e ({') is the solution to the 
following boundary value problem 



where as usual rj = 77(e). 

We finally fix positive constants 01, 02, q, /?, and C2, < Ci < C2 < 1. We 
also define the following space 1^(0) ={»£ W^ 1,2 (fi) : w = a.e. in fli}. To 
any w £ W(V£) we associate the function v = 1 — v. We remark that v £ M /1,2 (J7) 
and v — 1 almost everywhere in Oi . All these constants and the notation will 
be kept fixed throughout the paper. 

4 Determination of material losses 

In this section, the main of the paper, we shall consider the problem of deter- 
mining material losses. We begin by defining suitable classes of material losses. 

Definition 4.1. Let us fix a positive constant 5. We say that B is an admissible 
class of material losses if the following holds. First, any K £ B is an admissible 
defect such that dist(-R", > 5 and Gk is a domain with Lipschitz boundary. 
Second, we assume that, for some constant C, we have T-L n ^ 1 (8Gk) < C for 
any K £ B. Finally, we assume that the set {Gk '■ K £ B} is compact with 
respect to the Hausdorff distance. 

In the remaining part of this section, let us fix B, an admissible class of 
material losses in the sense of Definition 14.11 We assume that the unknown 
defect Kq belongs to B. We observe that, as in Proposition 12. 51 we have there 
exist a constant q > 2 and a constant C > 0, not depending on /q, such that 
Vu G L q (Q,R N ), in particular 



Here the constants q and C depend also on s and on Kq. In the sequel of the 
section, we shall fix q > 2 as such a constant, which depends on Kq, among 
other things. We define qi = (q—2)/(2q) and we observe that < qi < 1/2. We 
also define the following set. For any positive constant a, we say that v £ H(a) if 
v £ W 1 ' 2 ^!, [0, 1]), v = 1 almost everywhere in Cti and there exists K £ B such 
that v > C2 almost everywhere in £l\B a (Q\G k) and v < c\ almost everywhere 
in 0\B a (Gi<-). We observe that, by the compactness of the class B with respect 
to the Hausdorff distance, such a set H(a) is closed with respect to the weak 
W 1 ' 2 ^) convergence. 



(3.10) 



{ 



div(w v (v)'Vue) = in £1 
w r) (v)Vu s • v = / e on dQ, 



||Vwo|U<r(n) < C||/ ||l-(7). 
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For any e, < e < 1, we define Q £ : W(Q) — > R as follows. For any v G W(ft), 
recalling that v = 1 — v, we set 

(4.1) G s (v) = ^ j \u £ -g £ \ 2 + b [ «;„(«)|Vu e | 2 + i [ W(v)+r][ \Vv\ 2 . 



Here r\ = r](e), o v — o v (2), w v = w v ^(v) = tp v ( e )(v) and u £ = u £ (v) is the 
solution to (|3.10l) . Here and in the sequel of this section, we may also set the 
constant b — 0, that is we may drop the second term of the functional. 

Then, for any e, < e < 1, we define Q e : L 1 (il) — > [0, +oo] as follows. For 
any v € we set 

(4.2) g £ (d)=g e (v) if v € W(ty and v = (1 - v) € H(a e ), 

whereas G £ {v) — +oo otherwise. 

Theorem 4.2. Besides the previous notation and assumptions, let us further 
assume that the following constants satisfy < (3 < q < 2, and that 

r]{e) 2qi 
limsup ; — < +oo, 

and, finally, that a e > 277(e). 

Let uq = u{f(s,Ko). For any e, < e < 1, let 

m £ = mt{g e (v) : v € L\Q)}. 

Then we have that, for some constant C, m £ < C for any e, < e < 1. 
Furthermore, if v £ G L 1 (0) is such that 

Qe{v E ) < C for any e, < e < 1, 

the following holds. For any e, let v £ = l — v £ and u £ = u £ (v £ ). Then we have that 
ip v (e)( v e)u £ — > uq strongly in L p (fl) for any p, 1 < p < +00, and ip v ( £ ){v £ )Vu £ 
converges to Vuo weakly in L 2 (il). 

Furthermore, for any constant c, c\ < c < C2, the sets {v £ > c} converge, as 
e — > + , to Gk in the Hausdorff distance. 

Remark 4.3. We remark that the theorem in particular hold for a family 
v £ G L 1 ^) of minimizers or quasi-minimizers, that is satisfying 

lim (Q s (v e ) - m £ ) = 0. 

Proof. By Proposition 13.21 we infer that there exists a constant C such that 
for any e, < e < 1, and for any v € W(fl), we have 

w v {v)\\7 U£ (v)\ 2 < C and \\u e (v)\\ L oo (n) < C. 

By a construction pretty similar to the one used in Proposition 4.5 in |23j . 
we may construct v £ € W(fl) and u £ £ W 1 ' 2 ^), for any e, < e < 1, such that 
the following properties hold. For any e, < e < 1, first v e = (1 — w e ) G H(a £ ) 
and 

- / W> £ ) + 77 / |Vt; e | 2 < C. 
77 Jo 
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Second, on 7 we have u e | 7 = uo| 7 = go- Finally, 

,(iU|V( W£ -u e (5 £ ))| 2 <C^. 
By Poincare inequality in we conclude that 



<> 



(4.3) / \u £ (v E ) ~ 9eV < 2 / \Ue(Ve) ~ u s \* + / \g - g £ \ 2 < C(e« + e 2 ). 

We immediately conclude that for a constant C we have m £ < C for any e, 
< e < 1. 

For any neN, let e„ > be such that lim„ e n — and let Q n — Q En . Let 
v n be such that G n (v n ) < C, for any n € N. Then by Remark 12. 31 we obtain 
that, up to a subsequence, u n converges in L 1 (J7), and actually in L p (il) for any 
p, 1 < p < +00, and almost everywhere in fi, to a function u. Such a function 
i> is such that P(v) is finite. Furthermore, by the definition of H(a) and the 
compactness properties of B, we may also assume that there exists K G £> such 
that v = 1 almost everywhere in Gif and v = almost everywhere in f2\Gx- In 
other words, v = XGjc • 

Let us call w n = w n (e n ){^n) and u n = u £n (v n ). Let us notice that ^w n \/u n 
is uniformly bounded in L 2 (Q,M. N ), therefore, up to a subsequence, ^Jw n \/u n 
converges to V € L 2 (tt,R N ) weakly in L 2 {Q,R N ). 

Since v n — > xg k almost everywhere in f2, we conclude that also w n and ^/w n 
converge to xg k almost everywhere in VL and in L p (tt) for any p, 1 < p < +00. 

For any B2Ft{y) C {VL\Gk), we have that u>„ converges to zero almost every- 
where in B2n(y). By the uniform L°° bound on u n and by the dominated con- 
vergence theorem, we conclude that Jb 2R ( v ) w n^n -> as n 00. By the Cac- 
cioppoli inequality described in Lemma 12.71 (|2.14p . we conclude that ^w n V u n 
converges to strongly in L 2 (Bu(y),R N ), consequently V = almost every- 
where in £1\Gk- We conclude that w„Vh„ weakly converges to V in L 2 (f2, R w ) 
as well. On the other hand, again up to subsequences and by using the property 
of H(a), we may follow the arguments of the proof of Proposition 4.3 in [53] in 
order to find a function u with the following properties. First, y/w n u n converges 
to u almost everywhere in f2 and consequently in L p (f2) for any p, 1 < p < +00. 
Second, by the same reasoning above, we conclude that it = almost every- 
where in Q\Gk and that also w n u n converges to u almost everywhere in SI and 
in L p (Vl) for any p, 1 < p < +00. Then, we have that u € W 1,2 (Gk), u is 
harmonic in Gk and Vu — V in Gk- We also have that on 7 and on 7, u n 
converges to u strongly in L p (^f n 7) for any p, 1 < p < +00. As a consequence, 
u = g on 7. 

Let us take any function if £ W 1,2 (Gk)- Since Gk is a domain with Lip- 
schitz boundary, therefore it is an extension domain, we can find a function 
cp £ W 1,2 {VL) such that (p = f on Gk- We conclude that for any n £ N we have 



Since, as n — > 00, 



w n Vu n V<f= f e „<f. 
n J7 



id„Vu„ ■ V if —> / Vu • \/tf 
a JG K 
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and 

/ fe n <P ~> / fo<P, 

J 'y J ~y 

we conclude that 

/ Vu-Vp = [ f (p for any ip € W 1,2 (G K ). 

JG K Jl 

Then u solves p.l[) - (l3.2[) - p.3p with K replaced by K. Then, since u = u Q 
on 7, we conclude by using Theorems 3.3 and 3.6 in [21J that u — uq almost 
everywhere in f2 and that Gk = Gk - The rest of the proof easily follows. □ 

We conclude this section with the following existence results, which may be 
easily proved by the direct method. 

Proposition 4.4. The following problems admit a solution. 
(j) mmQ £ on W(Q), with constraint < v < 1. 

(ii) minC? £ on W(il), with constraints < V < 1 and v G H(a E ) {that is there 
exists the minimum of Q e over L l {£l)). 



5 The crack case 

In this section we shall deal with the determination of general defects, in par- 
ticular of cracks. We begin by recalling results proved in [33]. We include them 
here for the convenience of the reader and to compare them with the new results 
devoted to the determination of material losses, in particular of cavities, which 
we treated in Section For what concerns the classes of admissible defects we 
shall use in this section, let us begin with the following definition. We limit 
ourselves to the two or three-dimensional case, however it is not difficult to see 
how these definitions can be generalized to higher dimensions. 

If N — 2, fixed a positive constant L > 1, we say that T is an L-Lipschitz, 
or L-C 01 , arc if, up to a rigid transformation, T = {(x,y) € M 2 : —a/2 < 
x < a/2, y = tpi(x)}, where L _1 < a < L and ipi : M — > K is a Lipschitz map 
with Lipschitz constant bounded by L and such that ipi (0) = 0. For any a, 
< a < 1, we say that T is an L-C l a arc if ipi is C 1,a and its C 1,a norm is 
bounded by L. The points (a/2, tpi (a/2)) and (— a/2, ipi(— a/2)) will be called 
the vertices or endpoints of the arc T. 

Let us consider now the case N = 3. Let T be the closed equilateral triangle 
which is contained in the plane LT — {(x, y, z) e R 3 : z = 0} with vertices 
Vi = (0,1,0), V 2 = (-73/2,-1/2,0) and V 3 = (V3/2, -1/2, 0) and T' C M 2 
be its projection on the plane II. Fixed a positive constant L > 1, we call 
an L-Lipschitz, or L-C 0,1 , generalized triangle a set T such that, up to a rigid 
transformation, T = {(x, y,z) G M 3 : (x,y) € <p(T'), z — ipi(x,y)}, where 
: M 2 — > M 2 is a bi-Lipschitz function with constant L such that (f(0) = and 
ipi : M 2 — > M is a Lipschitz map with Lipschitz constant bounded by L and such 
that (fi(0) = 0. For any a, < a < 1, we say that T is an L-C 1,a generalized 
triangle if (pi is C 1,a and its C 1,a norm is bounded by L. 

In both cases, the image through <p of any vertex or side of T" will be called 
a generalized vertex or generalized side of tp(T'), respectively. The image on 
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the graph of ipi of one of the generalized vertices of <p(T') will be called a 
generalized vertex of F, whereas the image of one of the generalized sides of 
<p(T') will be called a generalized side of T. We also remark that there exists a 
constant L\ > 0, depending on L only, such that we can find </?2 : R 3 — > R 3 , a 
bi-Lipschitz function with constant L\, such that T = <p2(T). 

Definition 5.1. Let us assume that f2 C -Br C M. n , with R > 1 and A" = 2,3. 
For any positive constants L > 1, S and c, c < 1, any k = 0, 1 and a, < a < 1, 
such that fc + a > 1, we define B(N, (k, a), L, S, c) in the following way. We say 
that A £ B(N,(k,a),L,5,c) if and only if A C -B2J?,, there exists a positive 
integer rt, depending on A, such that A = 1J,™ =1 1\, I\ an L-C k a arc (if N = 2) 
or generalized triangle (if AT = 3) for any i — 1, . . . , n, such that the following 
conditions are satisfied: 

i) for any i,j £ {1, . . . , n} with i 7^ j, we have that either I\ n Tj is not empty 
or dist^i,^) > 5; 

ii) for any i, j £ {1, . . . , n} with i 7^ j, if I\ n Tj is not empty then I\ n is 
a common endpoint V if N = 2 and either a common generalized vertex 
V or a common generalized side 7 if A" = 3. Furthermore, in such a case, 
for any x £ T$ we have dist(a;,r :; ) > c\x — V\ or dist(x,rj) > cdist(x,7), 
respectively. 

Let us remark that there exists an integer M, depending on AT, R, L, 8 and 
c only, such that for any A £ B(N, (k, a), L, 5, c) we have that n < M. 

More importantly, we have that any of the classes B described in Defini- 
tion [5TT] is non-empty, is composed of non-empty compact sets and it is compact 
with respect to the Hausdorff distance. Finally, if A belongs to any of these 
classes, then / H N ^ 1 (A) is bounded by a constant depending on the class only. 

For the time being, let us fix B as one of the classes of Definition 15.11 We 
call the constant k, a, L, 6 and c the a priori data related to B. For any such 

class B we call B' the class of admissible defects K such that dist(A', fii) > 6, 
H N ~ 2 (KndQ) < +00 and there exists A £ B such that K c A and H N ~ 2 (Kn 
A\K) < +00. 

Moreover, we say that K £ B' satisfies Assumption A if the following holds. 

Assumption A. We assume that, for any xo € K fl O, there exists r > 0, 
depending on xq, such that for any U connected component of (Q\K) n B r (xo) 
we can find n > 0, an open set Ui, such that U n B ri (x Q ) C Ui C U, and a 
bijective map T : U\ — s- (—1,1)-" such that the following properties hold. The 
maps T and T -1 are locally Lipschitz and there exists a constant C such that 
||DT|| and ||DT _1 || are bounded by C almost everywhere. By the regularity of 
Q = (—1, 1) , T _1 can be actually extended up to the boundary and we have 
that T _1 : [—1, 1} N — > M is a Lipschitz map with Lipschitz constant bounded 
by C. Furthermore, if we set T = [-1, l]^ 1 x {1}, we require that T -1 (r) = 
dUi n K , T-^O,...^,!) = x and T-\y) £ Q\K for any y £ [-l,l] w \r. 

We assume that, for any x n £ K n dVL, there exists r > 0, depending on 
xq, such that for any U connected component of (Cl\K) n B r (xo) we can find 
r% > 0, an open set Ui, such that {/ n B ri (xo) C J7i C U, and a bijective 
map T : U\ — > (0, 1) x (—1, l)^ -1 such that the following properties hold. The 
maps T and T _1 are locally Lipschitz and there exists a constant C such that 
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\\DT\\ and ||Z?T _1 || are bounded by C almost everywhere. By the regularity of 
Qi = (0,1) x (— 1, l) w_1 , T -1 can be actually extended up to the boundary 
and we have that T _1 : Qi — > M. N is a Lipschitz map with Lipschitz constant 
bounded by C. Furthermore, if we set Ti = [0, 1] x [—1, l] N ~ 2 x {1} and T 2 = 
{0} x [-1, l] N ~ l , we require that T" 1 ^) = dU x n K ± T- 1 (T 2 ) = dUi H dVt, 
T- X (Q,...,Q,1) = x and T~ l (y) G fl\K for any y G Qi\(Ii UT 2 ). 

In the sequel we shall fix positive constants L > 1, S and c, c < 1, and a, 
< a < 1. We also assume that 51 C £?r, for some fixed constant R > 1. Let 
6 = $(jV, (1, a), L, 5, c). We assume that the unknown defect K$ belongs to B 1 
and that it satisfies Assumption A. We recall that examples of defects satisfying 
Assumption A are described in [22] [23] . 

The next proposition states that the gradient of uq satisfies a higher inte- 
grability property. 

Proposition 5.2. Under the previous assumptions, there exist a constant q > 2 
and a constant C > 0, which do not depend on fo, such that Vuo G £ g (f2, 
in particular 

||Vu |U*(n) < C||/o|U"(7)- 
PROOF. See the proof of Proposition 4.5 in [22]. □ 

We remark that the constants q and C in Proposition 15.21 depend also on s 
and on Kq. 

For any a > 0, we call -Hi (a) the set of functions v G W ,2 (&l, [0, 1]) such 
that v = 1 almost everywhere in fii and for some A G B we have w > c\ almost 
everywhere in £l\B a (A), where again c\ is a constant such that < c\ < 1. 

For any < e < 1 and any g > 2, let us define : W^' 9 (fi) x ->■ K 

as follows. For any (u,v) G W^' 9 (f2) x W(fl), recalling that v = 1 — v, we set 

(5.1) JFf (u, 5) = °±\u - u e \ 2 Wn + % j \u - g £ \ 2 + 

b [ ^ n (v)\Vu\ q + - [ V(v)+n [ \Vv\ 2 . 
Jn V Jn Jn 

Here r\ = 77 (e), o ,, = o r) (g), = u>,,( e )(fi) = t/tye)^) and u £ = u e (5) is the 
solution to (I3.10p . We also recall that 

\u-u e \l = [ ^ {£} (v)\V(u~u e )\ 2 = / V„( £ )(u)|Vu| 2 -2 / / e ii+ / / e « e . 

Then, for any < e < 1 and any q > 2, we define J-"* as the following 
functional on L 1 (f2) x L 1 (f2). For any (u,v) G i 1 (f2) x L 1 ^) we set 

(5.2) ^(u,5)=^(u,S) 

if (u,v) G W^ g (fi) x W(Sl) and u = (1 - v) G ffi(o e ), 

whereas F q {u,v) = +00 otherwise. 

Now we shall fix the constant q > 2 as the one defined in Proposition 15.21 
which depends on K , among other things. Again we set q\ = (q — 2)/(2q) and 
we observe that < q\ < 1/2. The following convergence result is the main 
result of [33]. 
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Theorem 5.3. Besides the previous notation and assumptions, let us further 
assume that the following constants satisfy 0<q<2,0<(3<2, and that 

limsup = — < +00, 

and, finally, that a £ > 2rj(e). 

Let uq = u(fo,Ko). Then there exists a constant Eq, Eq depending on s, fi, 
f2i, Q\, 7 and 7 only, such that for any E , Eq < E < +00, the following holds. 

For any < e < 1, let 

m s = in£{J*(u,v) : (u,v) £ L 1 ^) x L 1 ^) and \\u\\ L oo( U ) < E }- 

Then we have that, for some constant C , m e < C for any < e < 1. 

For any n £ N, let e n > be such that lim n e n = and let (u n ,v n ) £ 
L 1 (f2) x L 1 (f2) 6e suc/i f/ia£ ||un||i<»m) < -B and 

T1 n {u n , v n ) < C for any n £ N. 

Then, up to a subsequence, u n u strongly in L p (fl) for any p, 1 < p < +00, 
and ^(e„)("n)V« n — > Vu strongly in L P (Q) for any 2 < p < q, where u = uq 
almost everywhere in Gk and Vu = V«o almost everywhere in Q. 

Furthermore, there exist compact sets A C £1 and A £ B, such that A C A 
and T-L N ~ 1 (J {u)\A) — 0, satisfying the following property. For any constant c, 
< c < c\, the sets {v n < c} converge, as n — » 00 , to A in t/ie Hausdorff 
distance. 

An analogous to Proposition ^. 4l holds true, again easily proved by the direct 
method. 

Proposition 5.4. Let Eq be as in Theorem 15.31 Then for any p, 2 < p < q, 

and any E, Eq < E < +00, the following problems admit a solution. 

(i) min T% on W^ p (fl) x W(Q), with constraints < v < 1 and \\u\\loo(q) < 
E. 

(ii) min.Ff on W^' p (f2) x W(Q), with constraints < V < 1, v £ Hi{a e ) and 
\\ u \[L°°{n} < E [that is there exists the minimum of J- p over L 1 (fi) xL 1 (fi) 
with the same L°° bound on u). 

Let us now consider the main differences between the cracks and material 
losses cases. Our aim is to show the optimality of Theorem [531 by showing that 
a reduction to a functional depending on the phase-variable only, with similar 
convergence properties, may not be feasible. As we have shown in the previous 
section such a reduction is instead possible in the material loss case. 

By Proposition [53] and Proposition [321 we infer that there exists a constant 
C such that for any e, < e < 1, and for any v £ W(fl), we have 

w v {v)\\/u e (v)\ 2 < C and \\u E (v)\\ Lo , {n) < C. 

Furthermore, there exists q(e) > 2, depending on N, fl, s and e only, such that 
u e (v) belongs to W 1 ' 11 ^ {£1). We can also find a constant C\, depending on N, 
7i s > ll/o||,L s (dn) and e only, such that for any v £ W(tt) 
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We remark that the dependence of q(e) on e is through o v ^ and that, unfortu- 
nately, it might happen that q(e) — > 2 + and C\ — > +00 as e — > + . 

Let us consider the following operator. For any e, < e < 1, we define 
H £ : W(Cl) -> W^ 2 (n) as follows 

Heiv) — u E (v) for any v € W^(f2). 

We recall that for any r, 1 < r < +00, we endow M^' r (f2) with the norm 
ll u llw l r (n) = ll^ u lli r (fi) f° r an y u G W^ ,r (r2). We observe that T-L e is continuous 
with respect to the weak-W 1,2 (£l) convergence in W^fi) and strong convergence 
inW^ 2 (0). 

We obtain that for any g, 2 < q < q(e), we have that H E : W(£l) -!• 
and that for any q, 2 < q < q(e), Ti. e is continuous again with respect to the 
weak-W 1 ' 2 ^) convergence in W(f2) and strong convergence in W^' q (fl). 

Then for any q > 2, let us dchne F q : W(£l) — > [0, +00] as follows. For any 
v e W(ft) we set 

(5.3) fl(v)=F2(H E (v),v) = 

% f \u e (v)-g £ \ 2 + b f ^ n ( v )\Wu £ (v)\ q + - [ V(v)+n f \Wv\ 2 . 
e p J 7 Jn V Jn Jn 

Let us notice that for any q > 2, we have that there exists min on W(il) with 
the constraint < v < 1, and with the constraints < v < 1 and v € H\{a e ) as 
well. 

We investigate whether, for some q > 2, we may have convergence properties 
for .FJ as we have for FJ. We observe that Q e is equal to P 2 but to replace the 
single- well potential Vwith the double- well potential W. It would be desirable 
to have a convergence result for P 2 , or at least for P% with some q > 2, as we 
have for Q £ , Theorem 14.21 By counterexamples we show that difficulties arise in 
both cases. We begin with the case q = 2 and then we deal with the case q > 2. 

By the construction used in [231 Proposition 4.5] the next proposition im- 
mediately follows. 

Proposition 5.5. Under the assumptions of Theorem \5.3\ and ifO < j3 < q < 2, 

we can find v e for any e, < e < 1, such that the following holds. For any s, 
< £ < 1, we have, first, that 

f 2 e {vs) < C. 

Second, {v £ < 1/2} = {x € : dist(x, -Ko) < £ v + r//2} where £ v = y/r]o v . 
Finally, for any n € N, let e n > be such that lim n e„ = and let v n ~ v £n 
and u n = u en (v n ). Then, up to a subsequence, u n — > u strongly in L P (VL) for 
any p, 1 < p < +oo, and ip n ( £n ) (v n )\?u n — > Vu strongly in L 2 (Q), where u = uq 
almost everywhere in Gk cind V« = V«o almost everywhere in fi. 

In terms of T- convergence, we have obtained a kind of r-limsup inequality. 
What is missing is the corresponding r-liminf inequality, because taking q = 2 
does not guarantee enough compactness. In fact the solutions to the correspond- 
ing weighted elliptic problems may converge to a function which is not a solution 
to a material loss direct problem, as we shall show in Example 15 .61 where we use 
the instability of the Neumann problem with respect to boundary variations. 
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In any case, trying to solve the inverse problem by minimizing on W(fi) 
with the constraints < v < 1 and v € H\{a e ), might be a good strategy. We 
recall that in this case the assumption < (3 < q < 2 should be adopted. In 
fact, minimizing T 2 is numerically simpler than minimizing and still leads 
to good numerical reconstructions. In fact this method is adopted in [20] and 
the numerical simulations presented there show its efficacy. 

Example 5.6. Let us consider the following example. Let D be a smooth 
bounded domain of l^" 1 , N > 2, and let A 2 > be a Neumann eigenvalue 
for -A on fl and let / be a corresponding eigenfunction, that is 

-A/ = A 2 / in D 
V/ • v = on dD. 

We notice that J D f = and we may normalize / in such a way that J D |/| 2 = 1. 

For some constant T > 2, to be fixed later, let ft = D x (0,T) and let 
G Ko = D x (0, 2), that is K = D x {2}. Let 7 = 7 = D x {0}. Then let u be 
a solution to 

Aug = in Q 
Vw ■ v = / on Z) x {0} 
Vuo ■ f = on iio 
Vu o -f = ondDx(0,T). 

We normalize uo in such a way that J" Uq = and, by separation of variables, 
we have that 



1 ^ /fa 



A 



cosh(2A) . . . 



x € D, ye (0,2), 




whereas ito may be chosen identically equal to in D x (2,T). 

By a simple computation, again by separation of variables, we may find T > 2 
and > and two functions u~ and u + such that the following conditions hold. 
First, u~ = uq in D x (0, 1) and u + solves 

= infix (1,T) 
z/ = on D x {T} 
v = ondDx(l,T). 

Second, the following transmission condition holds true on D x {1} 

u~(x, 1) = Uy(x, 1) = fi(u + (x, 1) — u~(x, 1)), x E D. 

By following Q15], we may then construct a Neumann sieve C Kq, S > 0, 
such that if solves 

Aus = in f2\Jf 5 
Vu^ ■ v = f on 7 
Vuj • 1/ = on 9(0\if 5 )\7 

the following holds. We have that, as 5 — > + , ug converges to u~ = uo weakly 
in H 1 (D x (0, 1)), and strongly in L 2 (D x (0, 1)), and us converges to u + weakly 
in H X {D x (1,T)), and strongly in L 2 (D x (1,T)) . Therefore, the Cauchy data 
of ug on 7 converges, for instance in L 2 (^f), to the Cauchy data of uq on 7. 



21 



By using Proposition 15.51 to approximate K g and ug, for any n <E N we can 
find e n > 0, rj n > and v n such that 

Fe n (vn) < C for any n € N, 

and that, asn-> oo, the following holds. First, e n —> + and rj n — > + . Second, 
if w n = w eii (t>„), ueN, then we have that u n converges to u" — uq strongly in 
L 2 (D x (0, 1)) and u n converges to u + strongly in L 2 (D x (1,T)). Furthermore, 
Vu„ • z/| 7 = Vuo • i^| 7 for any n € N and 

1 1 fin — m o||l 2 (7) -^0 as n — >■ oo. 

Therefore, even if the Cauchy data of uq on 7 are well approximated by 
those of u £n (v n ), we have that v n = 1 — u„ is small in a region close to the 
corresponding Neumann sieve which is far away from the actual location of 
the looked- for defect Kq. This example shows also the difficulty in proving a 
convergence result without imposing any further condition on the region where 
v is small. 

On the other hand, one might try to minimize T% on W(fi) for some q > 2. 
If we take q > 2, then compactness and convergence would follow as a simple 
consequence of Theorem 15.31 but we may not guarantee that we can find a 
sequence of phase- field functions v n such that (v n ) is uniformly bounded. 

Again we use the constraints < v < 1 and v G Hi(a e ). If one would be able 
to find v e , < e < 1, such that J"|({i e ) < C for any < e < 1 for some constant 
C, then by Proposition 4.3 in [53], we would obtain the results of Theorem l5.3l 
replacing J"| with even allowing E to be equal to +00. 

We believe that constructing such functions v E for some q > 2 is a difficult 
task and that minimizing J 7 * for some q > 2 might lead to a not correct recon- 
struction. In Proposition 15.71 below we show the difficulty of obtaining such a 
uniform bound. 

In order to have higher integrability of the gradient of u E {v), we need to 
guarantee that w v (v) = ip n (v) is a weight satisfying certain properties, for in- 
stance those described by Stredulinsky in [3S]. An important class of weights 
for which these properties are satisfied is the so-called Muckenhoupt class A 2 . 

We recall that w, a non- negative measurable function over R^, is a weight 
if < w < +00 almost everywhere and w is locally integrable. We say that a 
weight w belongs to the Muckenhoupt class A2 if there exists a constant C such 
that for any ball B C M. N we have 




The best constant C for which (|5.4[) holds is usually referred to as the A 2 - 
constant of w. We observe that the ^-constant of w is always greater than 
or equal to 1. For more details about the Muckenhoupt weights and weighted 
elliptic equations, we refer for instance to [14] . 

Therefore, a reasonable assumption is to take w = Wn(v) belonging to the 
Muckenhoupt class A 2 and such that its A 2 -constant is bounded by C, for some 
fixed C. Without loss of generality we can assume that < w < 1 almost 
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everywhere and that w — 1 outside a given ball B2R. Consequently, we infer 
that there exists a constant C\, depending on C and R only, such that 



/ w- 1 < Ci. 



IB 

Proposition 5.7. Let us fix q > 2. Let e n , n £ N, be a sequence of positive 
numbers such that lim„ e n = 0. For any n € N, let = J 7 ^ and let v n G W(fl) 
be such that the following holds. For any n G N, we assume that < v n < 1 and 
we set n n = r)(e n ), v n = (1 - v n ) G Hi(a en ), w n = w n (v n ) and u n = u £n (v n ). 
For any n G N, we assume 



(5.5) / w~" < d 

and 



n 



W n \Wu n \ q + — I V(v n ) +7]„ f \Wv n \ 2 < C. 

n Vn Jn Jn 



Let us consider u as the solution to 
(5.6) 



Am = in Q 
Vu • v = /o on dQ. 

We assume that Kq satisfies the assumption of Theorem 15.31 and that u 7^ Uq, in 
particular that u| 7 ^ uq\j = go- We also assume, for the time being, that ft and 
/o are regular enough to guarantee that u G L°°(fl) and Vu G L°°(n,M. N ). We 
may also assume that actually v n provides a good approximation of Kq, namely 
that {v n < 1} C B an (Ko) where a n — a £n . 

Then we have that, as n 00, w n 'S/u n converges to Vu weakly in L (ft). 
Consequently, as n — > oo, we also have that J \u n — g £n \ 2 — > J \u — <?o| 2 ^ 
and J-^ n (v n ) — > +00. 

Proof. Let us compute 

w n \Vu n - Vu| 2 = / (f En - f )(u n ~u)+ / (1 - w n )Vu ■ V(u„ - u). 
n Jon Jn 

By the uniform bound on u n , we easily obtain that 

(fe n - /oH'Sn - asn^oo. 

/an 

Let us evaluate the other term. We have 

f (1 iv ) 
Jn -fu^n 

We apply Holder inequality with coefficients q, p and r such that q^ 1 + p^ 1 + 
r _1 = 1 and we obtain 



(1 — w„)Vu ■ V(u n — u) < 

l-w n \P\Vu\P) . 



CM '(X m " lv(i "- i,|a ) '(/„' 
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We use our assumptions to infer 



/ {l-w n )Vu-V(u n ~u)<c( ( w- r ^) 1 ( [ \l-w n \ p 
Jn \Jn J \Jn 

with C independent of n. We may choose r such that < r/q < 1, therefore, 
since < w n < 1, we have Wn T ^ q < w^ 1 . Hence, by (|5.5p . we conclude that 
J n w n \\7u n — Vu| 2 goes to zero as n — > oo. We obtain that / \u n — u\ 2 goes to 
zero as well. We then apply Theorem 4.4 in [23] and the proof is concluded. □ 

Therefore, even if v n is a good phase-field approximation of Kq, u n is not a 
good approximation of uq. In order to have that u n approximates u$, we need 
to require that v n is very small close to K , in such a way that violates (|5 . 5|) . 
In turn, this might suggest the fact that higher integrability and the correct 
approximation might in some sense oppose each other. 

Let us conclude by observing that (|5.5|) is a kind of minimal condition to have 
J n w n \S/u n \ q uniformly bounded. We wish to point out that potential theory for 
weights whose inverse is not integrable has been developed, see for instance [TO] 
for the case of weights w — Lo 1 ^' p / N , 1 < p < N, where w is a so-called strong 
^4oo-weight. Strong ^loo-weights have been introduced in [5]. Following [23], an 
important example of strong ^oo-weights is given by 

u(x) = min{l,dist(x, A) s }, xER N , 

where s > and A is a suitable compact set. In [23] Proposition 4.4] it is shown 
that w is a strong ^loo-weight for any s > provided A is uniformly disconnected. 
On the other hand, no strong Aoo-weight may vanish on a rectifiable curve, 
therefore this class of weights seems to be not apt to approximate hypersurfaccs 
as we require in our application. 




6 Differentiability of the functionals 

In this last section, we investigate the differentiability properties of J-j, Q e and 
JF%, for a fixed e, < e < 1, and any q > 2. For this purpose, we further assume 
that the functions ip, V and W are actually of class C 1 and such that their 
derivatives are bounded and uniformly continuous all over R. 

We define the following spaces. For any p, 2 < p < +oo, let us call L p (il) = 
{v £ L'P(n) : 5 = a.e. in £l x } and W p {Vt) = W 1,2 (ft) H L p (fl), with norm 
\\v\\L p (n) = \\v\\lp(Q) and ||w||w p (fi) = \\v\\Lp(n) + \\Vv\\l*(si)- To any v € L 2 (fl) 
we as usual associate the function v — 1 — v. If v belongs either to L P (Q,) or to 
W p (ft), then v E L p (il), v = 1 almost everywhere in f2i, and, provided < v < 1 
almost everywhere in f2, we also have < v < 1 almost everywhere in O. We 
observe that W^(^) = W(fl) as previously defined. We also recall that W^ q (£l) 
is equipped with the norm IMI^i.g^ = ||Vtt||i«(n) for any u € W^ ,q (Vl). 

We recall that for any e, < s < 1, we define W E : L 2 {£1) -> W^' 2 (il) as 
follows 

H E (v) = u e (v) for any v e L 2 (fl). 

It can be shown that for any vq € L2(£l) such an operator T-L e is differentiable 
in v with respect to the L°°(r2) norm. Let DH £ (vq) : L^ft) — > W}.' 2 (Vt) be 
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the differential in Vq. Then for any v in L^^) we have 

DU e (v )[v] = U £ (v Q ,v) 
where U e = U £ (vq,v) £ W^' 2 (tt) solves the following problem 

f dWty v (v )VU e ) = div(^(« )6V(?4(«o))) in SI, 
1 J \ ip v (v a )VU E -u = ondfl. 

Here, obviously, vq = 1 — 5q< 

We recall that for any vector valued function G £ £ 2 (f2, div(G) defines 
a functional on W 1 ' 2 (SY) in the following way 

div(G)[</>] = - I G V(f> for any € VF 1 ' 2 ^). 

Therefore, the weak formulation of (|6.1I) is looking for a function [/ E g W^ ,2 (f2) 
such that 

il> v (vo)VU e 'V<p= [ ^{v Q )vV{U e {vo))-V<p for any ip £ W^ 2 (Sl). 
Jn 

Here, and analogously in the sequel, the differentiability has to be understood 
in the following sense. For any v in L^iSl) 

U e {v +v) = U E {v ) + DH £ {i Q )[v] + R(v) 

where 

\\ R (v)\\w^(Ci) . 
hm — -7— = 0. 

||«IUo<.(ii)->o ||«||i~(n) 

For any q > 2, let us consider the functional f§ : W}j' q {Sl) x W(SL) -> K. For 
any (uq,Vo) £ W,J' 9 (f2) x W(O), fij is differentiable in (uo,5o), with respect to 
the W^{Sl) x WooCfi) norm. Let DJ|(w ,u ) : x W^ft) ^ M be the 

differential in (wo^o)- Then, for any (w, w) g W^ ,9 (f2) x Woo(^)j we have 

(6.2) DF2(u ,v )[(u,v)] = 

% / (2^(bo)Vuo ■ Vu-T//> )|Vii | 2 5) + 
£V Jn 

(g^(^o)|Vwo|^ 2 Vu • V« - ^Oo)|Vu | 9 u) + 
?7 Jn 

With the same computation, we infer that the functionals F 2 : W{Sl) — ► R 
and C7 e : PT(fl) — > R are differentiable in -Do for any vo £ W(Si), with respect to 
the Woo(fi) norm. Let DF 2 (v ) : W^tt) -> R and DG £ (v ) : Woo(O) ^ M be 



(-V / («o)u)+ 277 / V£ • ve. 
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the differentials in (uq,vo). Then, for any v £ W oo (0,), we have 

(6.3) Df^voM = J(n s (v )-9s)U s (v ,v)+ 

b f (2iPr,(vo)VH s (v ) ■VU s (v ,v) - ^(v )\VH £ (i )\ 2 i) + 
Jn 

- I {-V'(v )v) + 2ti I V5o-Vw 
V Jn Jn 

and 

2 f 

(6.4) DQ e {v )[v] = ~f {H s (v )-g £ )U e (v ,v)+ 

b I (2iP v (v )VH £ (v ) -VU e (v ,v) - V;M|VH e (£ )| 2 £) + 
Jn 

- [ (-W'(v Q )v) + 2 V [ Vw • VS. 
V Jn Jn 

It might be useful to have differentiability properties with respect to the 
W p (ft) norm, with p finite. In fact in this case W p (f2) is a strictly convex real 
reflexive Banach space and this is useful when we need to apply a gradient 
method in a numerical implementation, see |20] for details on the use of this 
information. In order to obtain such differentiability, let us now assume that ip' , 
V and W are Holder continuous for some exponent a, < a < 1, all over R. 

We recall again that, by Proposition 12.51 and Proposition 13. 2[ there exists a 
constant C such that for any e, < e < 1, and for any v £ Li2(£l), we have 

w v (v)\\7u e (v)\ 2 < C and \\u E (v)\\ L ^ {n) < C. 

Furthermore, there exists q{e) > 2, depending on N, Q, s and e only, such that 
u e (v) belongs to W 1 ' q< - £ \Q). We can also find a constant Ci, depending on N, 
ft, 7, s, ||/o||i»(an) an d s only, such that for any v £ L 2 (f2) 

\\Vu e (v)\\ Lq (e)(n) < Ci. 

We remark that the dependence of q(e) on e is through o n ^ and that, unfortu- 
nately, it might happen that q(e) — > 2 + and Ci — > +oo as e — > + . 

However, we may conclude that U e : L 2 (Cl) -> W^' s(e) (0) and its image is 
bounded in (Q). Furthermore, again by Proposition [231 we infer that for 

any do G £2(0) we may define as before DH £ (iio) and prove that DH £ (vo) : 
Lq(e)(q(e)+2)/(q(s)-2){ft) — * W-y 2 (O) is a bounded linear operator. 

Letp(e) = ■ Then, straightforward but lengthy computations allow 

us to show that for any vo £ L-2(£l), H e is differentiable in vo with respect to 
the LP( £ 1(n) and the W^ 2 (fl) norms. The differential is still given by (|63]l . We 
immediately infer that for any p, p > p{e), Q e is differentiable in vq, for any 
vq £ W(f2), with respect to the W P (SY) norm, with the differential given by 

USD. 

By an interpolation inequality, we may find qi(e), 2 < 171(e) < (q(s) + 2)/2, 
depending on q(e) and a only, such that for any q, 2 < g < qi(s), and any 
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Pi P > p( £ )i we have that, for any vq <E Z^C^), "H e is diffcrcntiablc in vq with 
respect to the L p (Vl) and the W^' 9 (f2) norms. Obviously the differential is still 
given by (|6.1[) . 

We conclude that for such q and p, and any 5 6 W'(Q), we have that J"| is 
is diffcrcntiablc in «o, for any -Do S W(il), with respect to the W p (Vl) norm. Its 
differential is given by the following formula. For any v € W P (Q) we have 

(6.5) D^(vo)[v] = / (H e (v ) ~ g E )U e (v ,v)+ 
Jn 

- [ (-V'(v )v)+ 2t) [ V^o-VS. 
V Jn. Jn 

An important final remark is the following. If N = 2, then we may actually 
choose p(e) — 2, and we observe that W2(f2) is a Hilbert space, with the scalar 
product J„ Vui • V-D2 for any Ci, U2 € If iV > 2, then it might happen 

that p(e) > 2 and therefore W p ( e ) (f2) has not a Hilbert space structure anymore. 
However, since p(e) is finite, W p ( e )(f2) is still a strictly convex real reflexive 
Banach space. 
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